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ADSTTACT 


10  Influence  of  :    it Ion  heat  transfer  on  the  fore 
convection  dominated  heat  transfer  problem  has  not  been 
established*  An  anal;  tic  solution  for  this  prob]   is 
presented  :?or  the  cc     ■■•'   radiation  an.1  convection  heat 
transfer     a  turbulently  flowing  hot  fir.*"  In  a  lor 
circular  pipe.  A  sample  numerical  calculation  Is  carried 
to  completion  for  pure  carbon  dioxide  ;;;a3  flowj    'n   a  t 
inch  diameter  pipe .  Constant  mean  value  fluid  properties 
arc  assumed  throughout  the  analysis*  The  rc3u?.ts  of  this 
restricted  solution  show  that  radiation  effects  do  cause 
a  ;aln  in  the  total  not  heat  transfer  rate  to  the  pipe  wall 
compared  with  an  equivalent  scro  radiation  solution.  ;!ov;- 
ever,  it  is  also  shown   that  the  sum  of  independent  radiation 
and  convection  solutions  cannot  be  used  for  an  adequate 
description  of  the  combined  radiation  and  convection  heat 
transfer  solution. 
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he  need  Tor  an  analytic  solution  to  this  combined 
radiation  and  convection  heat  transfer  problem  was  first 
brought  to  ray  attention  by  Professor  TVarron  ♦  "  ohsenow 
of  the  rtechonical  epartnent  of  "'•  .  .   is 

suggestions  and  counsel  led  mo  to  the  General  extension 
of  the  xTon     in  convection  heat  transfer  analy     or 
flow  to  include  t     Vects  of  radiation. 
Special  thanks  are  due  to  "rofessor  Koyt  C.  "lot to  1 
of  the  Chemical    Lneering  T*>partment  at  KIT  for  provid- 

the  method  of  so  lutlon  for  the  total  not  radiation 
heat  flux  received,  by  any  infinitesimal  control  volume  of 
gas  In  the  pipe*  ITis  consultation  on  many  of  the  other 
radiation  considerations  involved  in  th£a  thesis  was  also, 
of  special  value. 

dward  5.  Cohen,  graduate  student  in  Chemical  Engineer- 
ing at  fIT,  was  of  assistance  on  frequent  occasions  in  jiving 
me  instruction  and  advice  on  many  of  the  radiation  aspects 
of  this  thesis. 
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Heat  fcransfe:       *oble: 

-.per  at  ur  err-  have  become  more  'leant  recently 

^r  of  high  ter.porature  fluid  a  itions. 

her  t  ,  all  flow  heat  tra  •  - 

for  ^roblermc  involvi  *     the  hetcropol  ases  may 

ly  affected  or  control!  .tioi  .  :;t 

noteworthy  gases  that  fit   into  this  category  are  or 

and  carbon  ic!c« 

-o  analytical  i  '  an  ;  h  blea 
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anal        "      the 
review  and  study  of  the  independent  forced  convection  heat 
transfer  analysis  for  turbulent  floi  3  tubes, 

da  analysis  was  then  roooat  '     letail  with  the 
ion   of  t    :cor.sary  radiation  terms.  of 

this  solution  employed  the  analo    between  radial  tr;    jr 
of  heat  and  momentum  in  lor^;  tubes,  ae  carried  to  solution 

Is.  [l0| 
s  followin;       try  condition*     limitations  v?cro 
by  this  solution: 
1*   lamination  is  restricted  to  a  ,-ivoa  cross-seotion 
of  an  axisymetric  fluid  flow  in  a  Ion;;  circular 
tube, 
2,  "oderate  temperature  o::tre ~es  exist  so 

(a)  sical      *ties  •.a.,  bo  a    sd  to 
bo  independent  of  temperatures 

(b)  temperature  gradients  have  no  effect 
on  the  velocity  radientn. 

o  fluid  is  in  a  state  of  hstea    fully 
developed  turbulont  flow.  T:\o   velocity 

Lstribution  is  given  by  the  experimental 
data  of  N'ikuradse.  foj   page  154 
Tho  details  of  the  fir    fcase  (Phase  A)  of  the  analysis 
are  shown  in  Appendix  I.  The  steps  of  this  analysis  aro  out- 
lino*:'1  as  follows: 
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Stop  1 

cloct  a  differential  element  control  volume  a1 
cross-section,  vantage  of  this  axi 

element  sel      at  thia  stare  is  of  % 

<?r  by  ex  ac  shown      ures  "  an 

Step  2 

Carry  out  the  momentum  transfer  anal    i  in  order  to 
obtain  the  express:7    o^  the  apparent  shear  stress  at  a 
rncius.  The  momentum  equation  may  bo  sir        to  the 
following  '"' 
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Figure  II  a  -  Force3  actinc  on  fluid  within  control  volume. 
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Figure  II  b  -  Fluxes  of  enthalpy  and  kinetic  energy 
through  the  oontrol  surface. 
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Figure  II  c  -  Rates  of  shearing  work  at  the  control  surface 
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order  terras  and  contributions  of  relatively 

to  ob1    tho  foil         ■   equation  as  a  result. 
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►!•©  (  dq  )  is  the  not  radiation  heat  transfer  rate  to 
7*£  r 


this  control  volurvo  per  unit  volume  of  fas. 


Step  4 


Define  the  eddy  heat  transfer  diffuslvity,  &_,  tat 

v  T     5  -  £-JT  and  define  apparent  hoat  transfer  flux,  (q/A) 
r       ■         h*rj  app 

such  that  at  an;;  ■■  ;iven  radius* 


(q/A)aPp  5  (q/A)molocular  +  <<*/A>cddy    +   ^A>P 

(Forced  convection) 


(7) 
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(q/A)    a 
^  'app 


Co) 


where  (q/A)  is  a  function  representing  the  1  at 

transfer  rate  at  this  point  per  unit  area  normal  00   the  ax 
(This  term  is  evaluated  in  Phase  B) 
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•'.V      of  the  it  transfer 

IX |  (  q/'  )flt)T,j  HS«   '  '•  ^01  "C 

juirod  .     or  to  roli 

(a)  le  heat  tra  ifer       iienl;  h  I  (10) 

is  very  nearly  uniform  bh« 

(Padiatlon  heat  transfer  is  a  .^ant 

a  c-!iall  fraction  of  the  total  heat  transfer 
ill* ) 

(b)  For  a  given  problem  at  specified        .'  the 

ynolds  n     t     ra  i  -Itl  io 

ra<  rti    •  the  fluid j  to 

3  »ature  distributio     ch  is 


r of ore 


h  '  ■  0  (11) 


— r-   rr 

a  -v      ,  -0 

>r  the  case  c     form  wall  .     transfer  flux, 
(  q/A  )Qf  alo:i:  the  le     of  the  pipe;  this 
reduces  to 

(c)  The  time  average  valuo  of  the  axial  velo 
component j  vx  ,  is  assumed  to 
constant  and  Independent  of  r    i  ,  (slur;  flovr). 
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reasonably  correct  within  the  '  at  ions  of 
as  Lc  assit.      .;•• 


£to£  6 


Introduce  the  analogy  between  friction  an   ieat  transfer. 
Ton  the  momentum  aquations: 
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but  unfortunately  a  full  aT;alo,  :y  In  now  V  tho  ad- 
ditional radiation  term  In  the  heat  tran  equal;!  ..    is 
term  will  receive  more  cietailed  o  a.  [nation  during   iase 
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Those  expressions  are  combined  with  equations  (13)  a 
(14)  to  obtain  the  following  diffe    Lai  equation  relat 
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Intocrato  this  temperature  diet]       i  equation  to  find 


.  The  first  half  of  ■ 


the  temperature  t  as  a  function  o: 

this  equation  can  bo  integrated  directly  after  cert-. 

simplifications  aro        shown  in  the  v<    arman  . 

b   second  part  involving  the  radiation  heat  flux  will  require 
•r  :  rlc  integration,  and  '.7111  bo  discuss         o  • 

Integration  is  carried  out  separate"  the 

three  zones  which  are  used  to  describe  the  velocity   istri- 
butlon. 
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This  is  as  far  as  the  analysis  cai  bo  carried  without 
evaluation  of  the  radiation  heat  transfer  function  (q/A)p. 


In  order  to  evaluate  the  radiation  heat  transfer  cross.' 
tho  boundaries  of  the  given  rin^-shapec!  control  voluro,  it  is 
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convenient  to  focus  attention  on  .11 

cube  control  volume  1       -•   .    .is  cube  is  located  at 
d  cylindrical  coordinates  r^, 

ens ions  dr.  by  r-dG    r  dx- •    ~  let  radiation  heat 

-».        i. 

transfer  transmitted  bo  this  cube  may  be  divided  into  three 
components:   (1)   emission  from  the  cube,  {2)   absorption  of 
'iation  from  the  pi  >c   v/a;ls  and  (3)  absorption  of  radiatio 
from  al]     other  g  c. 

Unfortunately,  the  nature  of  the  radiatlor  pr 

it  reduction  to  a  relati\        e  two-  '    lonal 
on  as  in  the   case  of  the  convection  eval     - .   ?he 
and  down  stream  contributions  to  the  overall  radiate 
transfer  to  a  given  cube  are  significant  and  necessary  for  • 
adequate  evaluation. 

i3t,  the  absorption  of  radiation  from  all  the  other  ,- 
in  the  pipe  will  bo  considered.      ;as  is  assumed  to  be 
for  the  time  being  for  ,    " r:  -'.11  be  corrected 

later  on.   The  pipe  wall  is  assumed  to  bo  black. 

otep  1 

valuato  the  radiation  omission  (c~q),..  which  is  recoiv 
by  the  riven  cube  (1)  at  r  ,  -^  and  x-  from  Bone  other  cube 
(2)  at  distance  (z)  located  at  r0,  ftn  and  Xr>  as  shown  in 
ure  ITI#  The  description  of  the  evaluation  of  t 

:-ven  1        "or  the  general  case  by   .   .  Cohen,  [l] 
.  one  way  radiation  (d-q)r  from  cube  (2),  (dv0),  to  cubo  (1), 

-  ),  is    lal  to  the  o-  iss     »oru  cube  (2)  t    \ 
x. 


multiplied  1-v  the  tran.     nice  of  the 

(e    ),  and  by  tho  absorptl     of  cube  (1  , 

element  of  differ  ait  J  3  . 

(d4o)  =  C4KeEgdTg)  (KftdVx)  e  * 

•   ■' ' ■■  ■  —  '  *"  * 

4?Tb  (25a) 

Then  for  this  case: 


i         a  '   A 

where  E0  2  <g(T0)   •     $L  is  tho  emission  coefficient  :°or  tho  rag 

and  KL   is  tho  absorption  coefficient.  *  It   is  convenient  to 
rewrite  this   equation  in  d linens lonl ess  form. 

Let     a~«     CxsAr     (*oA)     -  2(r,A)(r*A)  cos- 
X  I    (xg/  ), 
and  dv,   a  r-,  drjdcUdx^ 
Then 
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Integrate  this  expression  along  the  complete  length  of 
the  pipe  In  order  to  obtain  the  radiation  c    Lon  (d  q) 
which  is  received  by  cube  (1)  from  a  thin  gas  r       it 
local.     r„  aic"  9  oxtc         full  len  >ipe« 

his  point  it  is  necessary  to  make  some  a  i  on 

teeming  the  axial  \  rature  variation  at  r  .     is 

jd   for  this  problem  that  t-ho  axial  t    rature  radient 
relative!;   tall  (and  constant  *  iv« 

cross-section  as  shown  in  'haae  ;.).  Therefore*  at  anj  ;iven 
ii  ~i    .  ,ion  the  temperature  at  any  given  be 

j  .od  cc  integrate  m. 

.  jroxo: 


^3 


a. 
71 


<fiX 


(27) 


It  is  i.ioro  Gonvc.xion^  ror  numerical  oval. 

s     SO 

le*        ^  s     coa^,  sin  ^  :  x 


v^ 


~~zr~i 


2  .      r 

f 


g 


-i 


r  Wr's 


a   •+■ 


Then 


u  cJ 


»^hj;)(£/X^<S(^^ 


K' 


(28) 
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This  Integral  can  bo  carried  out  numerically  for  a  suitable 

-  n  a 
ranr.e  of  valuos  of  (K.Ra),      plot  of  (o  duti(f  )  versus 

ciL 

at   constant  (K-Ra)  value  is  shown  in  Figure    .  The  Integrals 

of  these  curves  nultiplled  by  (  1    )  are  plotted  against  (a) 

a 

at  constant  values  of  (Kaft)  in  Figure  v.  These  curves  are 


then  cro3s-plotted  on  semi»log  paper, j 


3, 


c<_ 


versus  (  LR)  at  constant  (a)  in  order  to  aid  in  interpolation 

to  any  desired  value  of  {    }.   It  should  bo  noted  that  those 

a 

curves  can  be  used  for  a  wide  rang©  of  pipe  flow  radiation 
pre   .  s  and  are  not     ricted  to     specific     rical 
example* 

i?.ten  5 

o 

Determine  the  radiation  emission  (d**q)..  which  is  received 
by  cube  (1)  from  a  thin  cyl     .cal  shell  of         ' L  •  ■  the 
full  length  of  the  pipe  symetrically  located  about  the  center 
at  radius  r0.  The  previous  expression  is  them  to  be  integrated 
about  9_  at  constant  rs. 


CrfV3 


i^mfx^K^fj^)?  ^h\ 


(SO) 


e,;a 


I  Tote  that  (a)  is  a  know    iction  of  9p  for  any  given 


value  of  (  r-,  )  avid  (    )  .     is  therefore  most  co    Lent 

to  carry  out  this  Integral  at  sonc  specified  value  or  scries 
of  values  of  the  parameter  L  . 


;ures  VI,  VII,  and  VIII 
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Figure    IV 
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Figure  VI 


*s 
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Figure  VII 
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/*° 


Fipure  VIII 


a./ 
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'  3  accomplish  this    '  i,   ecu 

orocet 


(a)  plot 


i. 


CO 


J  <8^ 


versus  (q)  at  const a 


values       (Figure  V). 
(b)  c    '~.to  (a)  as  a  function  of  -r,   at  ;ivon  valuoo 


(c)   >lot 


versus  b0   a 


of  (  ^  }  and  (  ^2    )    . 

r   ^  ^ 

from  the  previous  plot  at  the  correspo  11      oes 
of  (a),  for  given  values  of  KJR,  _1  and  ;2  ♦ 


)  carry  out  numerical  integrations  of  the  result i 
plots  from  9 g  =  0  to  oo  -  Tf. 

The  results  of  the  I    ration  for  a  sanple  value 

— »  j»< 

r  .7°  are  shown  in   t  ure  v  plotted  anainst  (    )  at 

a  —  "■ 

constant  values  of  (  rl  ) . 


Step  4 

determine  the  radiation  (dq) ,  which    'eceiv< 
cube  (1)  froia  the  whole  gas  in  the  pipe.  That  is,  integrate 

the  proceeding  result  across  the  radius  (££  }  fron  0  to  1. 

r  ■  *  '\ 


K{f 


f?>0 


SLs  o 


ee  Figure  IX  for  a  sanple 


14   - 


Fl,,ur 
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Fl jure  X 


s.o 
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z  &■ 


' 


r 
ately,  the     Lstributioi     i  {  _£  ) 

la  usually  an  u      "or  any  given  proble  .      ■  111  le 

to  a  trial  and        tterati    »oeedure  to  e:      h  the 

proper  to  .  . 

p.  too  5 

icider  t.  :  of  variation  fro  :  the  : 

pipe  wall*  Plrstg  evaluate  the  radiation  emission  (d5q) 

Loh  Is  received  by  the  r;lven  cube  (1)  from  an  area  olo  x 

on  the  pine  wall  of       .ons  r-  6  >     ! 

o     o  o 

distance    (a   )  at  r     ■     ,   -G  a         3.    (The  pioo  wall   Is 

assuiTiod  to  bo  black* ) 

The  one  way  radiation,    (d'-'q).,   fron  the  wall  ar 

element  ^{dA   )   to  the  gas  cube    (1),    (dv, ^  is   equal  to  the 

gsion  fron  the  wa  .  iA   )•     ru     L?llo<     s     the 

>   w  o  ' 

cosine  o:    3  an^l :  ie  connecting         bhe 

nomal  to  dA  7(cos/9),  by  the  transnittance  of  the  inter- 
vening gas  (e  a  '' )  a  &  .~    the  absorptivity  of  the  ;:as 


cube   (1),     "V^l      . 
~   "o 

(r? )   ^ l JOfr~/*rw  e^Kf!l 

For  this  case,    (cosine  .3)   can  be  zlioxm  to   be    "" 


(r.ia) 


yTS  ft*-  ^rKcas9a  +  / 


*  See  Appendix  II 
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Thorofi 


(31b) 


Lch  can  bo  slmplifj  lenaionlees  fc  •>  l 


for  the  gas. 


I 


^^^Xt/XieXf') 


0- 


(l-ft!  c-sQe 


^fv<V 


So 


L 


Integrate  this  expression  along  the  complete  bh 

of  the  pipe  In  orde:  >btaIo  the  radiation  i   (d2< 

-  Is  rcjcivc^  bj  (1)  .  ;  vail  eleraont 

the  lengt:  .  ?.  -  ■     ■     ■  - 

o  o  * 

the  axial  variation  in  wall  tg  _,r 

:J03   of  this    integration. 


3:=°° 


;)j[^A,/);-Xt'X'-^-^W 


K  J  ^  ~* 


X  y> 


*s 


<y/xc 


Again  variables  are  changed  In  order  to  avoir'  infini 

Lot   o      & 

<  i  li •  +  1  -  2  rl  cos  a 

u     +  ■ "  >?  "         o 


TT^ 


am1. 


"<<■' 


=  a. 


o 


a   X 
o   o 


\  *  •// 


uf 


is   in'         il  car:  be  carried  out  ir  sally  for  "."oblc 

ra  f  values  of   (:L;,u    ),  lot  of      -*«&?  ,  . 

versus    ^   at  constant    (IC„Ra.)  values   is  shown  i 

ti       o 

9  interrals  of  those  curves  nxulti  by    £p     aro        Dtted 

ao 

a_]   at  com  ■  .  •  .         :oso 

o  u 

■   arc   ala  u.  -  . 

r      ?l* 

Loo    ^x      ^  '    *v*«-       _  oJ  j} ^'    versus    (ILJ  )   at  cc         tnt   (a   )• 


::,.  P  7 


ter  the  rad     bio  is  rece 

jube    (1)  from  the  whole  pipe  wall  b     3 

exprossio.         sund  the  arnfore  jog   i  o. 

a  •  7r 


«b 


1 


^E/)(^ 


K' 


r 


-*i* 


0-r^os<4i\  -^5 


J 


C"*  *        -O^      C*^' 


<£©     ( 


3  ■-  o 


A 


that   ■;,  s<^T~   is  constant  for  this   Inte/ratio 
is   outside  the    r   te^ral.      clowover,    (a   )   Is  .  Ion 


of  a 


o 


»  See  Figures         .  ,  ,   and  XV  and 
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Figure  XI 
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Figure 
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Figure   XIII 
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Fif-ure  XIV 
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Fi,oire  XV 
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Figure  XVI 
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Thoref ore,  the  same  procedure  an  used  In  ■:         ;11 


integration  shou!    \   earrlod  out 


(a)     plot 


A.  1 


^(- 


co  5 


<£<*¥. 


versu  at 

o 


constant  is  or  K  ...       ..} 

9 

(b)     evaluate  a^  as  a  function  of  -     at  given  values  of 
o 


r? 


(c)    iaot  (x  •  fi  coo  o 


"/j 


_L 


ISa-o 


t4    e.      _     *  CosV<tfy 


vorsui    - 


for  .^ivon  values  o:  and   (     1   ) 


•» 


a 


(d)  carry  out  numerical  integration  of     result in 

plots  from  *?  z   0  to  9„  "TT* 
o  o 

The  results  of  3  int(        fcions  at  as  of 

K  R  s  0.7G  arc  sho^jn  i         .  urc  slotted  (   rl    )• 

a  — 


Step  G 

,   determine  the  r.                             Ion  -(dq) 

the  cu'jo   {!)•     rhis  one  is  simple • 

(dq^     ft  *  4  K^^dr^^^j                                                 (36) 

or  in  dimension! ess 

(dq)      :  -  4(K  R)(K,R2)(dVl      )  7) 

€  O           1 


Stop  9 

Combine  those  radiation  components  from  e<    .ons  {LD), 
(55),  a     M  to   '  raine  tho  net  radiation  heat  tranaf 
tho  cube  (1) 


(dq)r  Z   (dq),.+  (dq)w4(dq| 


'ifture  XVII 
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Figure  XVII 
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Figure  XVIII 


J.* 


/.o 


s.o 


i  tonne  ahoul     be  corrected 

r^oai     lot  c.  i:rav  ar« 

should  be  multiplied  by  the  real  gas  trans...    se 
Tact  or  (x)  •  Thia  factor  may  bo  the     or  as  t 
of  the  total  on 

of  the  real  ../.b,  .  pa    i  (1  -  (x) 

represents  trans parent  no-  l-ab sorbin;  and  non-ai  Lttlng  ;as. 
(x)  is  stralut  .ion  oi^  t     I  at 

a  ;;ivon  ?  L  valuo  (  e  J  and  at  a  double  P.    ..Vac  (£  . 

2 

6 

!<=  •  € 

jrcfore,   for  tlio  real   -:as 

(dq)        =    (x)J  (x).(dqL      (::)    'do)  (39) 

t  c  e 

v  tha  total  net  radiation  heat  trans for  cross ins  the 
ndaries  of  this  cube  per  unit  of  its  area  nonaal  t 

P1^      —  /    d<*     \  z  /<lci  \   •   Since  t    jubi 
\r  d-       J 
V  1   1      /r 

tjpical  of  all  cubes  in  a      control  voluric  at  e 

radius  r.  ,  this  also  ::;    ;onts  the  net       on 

1 

heat  flux  to  the  rin   >ntrol  volume. 

This  now  represents  the  net  radiation  heat  flux  to  a 
jrol  volume  o  as  of  thickness  cl  {  ^1  )•  To 

obtain  the  not  rad5.ation  heat  flux  per  unit  dims 
trol  volume  thic     ,  simply      e  by  d  (  rl  )  • 
ves  the  roquirec  radiation  point  function,  (q/A)r, 


the  net  radiation  heat  transfer  to  a  control  volir.o  ;>or  unit 
area  normal   bo   the  ;/all  for  lslonless  u  losa  of 

control  volume* 

(q/A)     S      (q/A)  (q/A)      ,,-f    (q/A)      .      .  0) 

r  ras  rail  emission 


where  r* .  \ 


(q/A)       : 


[^x^k4e^i  fc  K*^>e=U^  <*> 


(qA)       i 

wall 


*"<  (^1(1  -*  «*J  f,\e  -^co,  v. *  &1  *     (4a, 


&-  :  c; 


(q/A)     ,     ,    s    -  4{x)r.(::  }  )  (45) 

emission  ©         1 


HIASE  C  -    furaerical  Temperature  Calculation  Tor  a  Lfic 

rbcample  '  " 


In  or£er  to  -1  bhe  application  of  tho  prace 

results  to  an  actual  problem*   a  specific  oxa  1  be 

solved  numerically . 

The  statement  of  the  -tto   3        is  as  foil 
First,   consider  a  non-radiatj  tical  7^- 

turbulently  flowing  in  o  lony  cool  pipe*  b  non-radiati 


a 

he  si  nificance  of  the  selection  of  this  terra  a.3   the 

radiation  function,  (q/A)_,  which  forms  a  part  of  (q/A)a<0D 
in  equation  (7)  is  discussed  in  detail  in  the  pp 
Discussion  of  Results. 
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an  di< 
co    tlon      b         letical  nan  to  the  tube  r:all  will 
bo  calculated  under  the  following  physical  conditions. 

(a)  Pip©  <  latter,  '    *  2  inches 

(b)  Reynolds  Number  of  flow, 

N   *  /?  VT>   «  20,000 
He    *%  ^ 
t  o 

(c;  Oas  pressure  p  ■  1  atmosphere 

(d)  Pipe  wall  temperature f 

T  «  550°H 
o 

(e)  Bulk  ras  tenpjrature, 

T.  ■  2000°^ 
b 

where  T,    rsiay  be  defined  as 

r 

-   •       [^    (V--VJr 

J.       ■     j   v  ''•»■' /^x.\  /■■  -  '  t/ [5]  (44) 


o 


c 


To  "  Tb        (    ^Cy.       ^  Jr 

a       'A'ln*. 

(f)  Oas  film  tenp 

All  physical  properties  are  evaluated  at  the 
film  temperature*  These  arc  tabulated  in 
Appendix  III. 

Prom  tbt  '■-  '  te  coefficient  of 

transfer  (h)  by  convection  only  is  determined!  fron  the 
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ant  on  number. 


ot 


V 


t  FT 
p  n 


i*.  • '   *- (45) 


/  +  J^U  [  6-  (-<■  %'  'V  SU)  +  fc^-*JUL] 


Then  the  net  heat  fl       j  wall  le  determined 

(q/A)   ■  h  (T   -  2.  )  (4G) 

O  O      D 

For  this  caso 

(q/A)   «  -  14,500  Btu/br,ftS 
o 

Next,  the  temperature  distribution  across  the  pipe  is  deter- 
mined by  use  of  the  von  Karman  method  of  solutio  . 
solution  is  identical  with  the  c     »d  heat  transfer 
temperature  distribution  solution  presented  by  equations 

(20),  (22)  and  (24)  when  (q/A)   =  0. 

r 

The  resulting  temperature  distribution  from  these 
equations  for  zero  radiation  is  shown  In  Figure  XIX, 

jw  shift  to  the  case  of  the  real  radiating  ;ao,  pure 
carbon  dioxide,  with  otherwise  identic.  9  properties. 
Assume  that  the  following-  physical  conditions  exist* 

(a)  Pipe  diameter  -  same  value,   ■  2   inches 

(b)  iber  of  flow  -  same  value, 

*  20,000 
Re 

(c)  Gas  pressure  -a         ,  ••■  ■  1  at 

( d )  ...         value    ■  550  F 
(o)  Total  •:        flux  for  t    as  at       II 

(now  Including  radiation)  -  sa      le,  (q/A)  * 

3. 
-  i '1,500  Btu/hr.ft   (not  the  same  as  the  nob 

heat  flux  to  the  wall  as  shown  in  I 
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Figure  XDC 


-r^i 


»•#• 


7tf>0 


/#«• 


*»<> 
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1 lUJ  *  .x.iZ 

te  '  .1  nor  11;  ,  so  :  '.'tios 

ie  bu] 

■  '.  '    .  • 

b 

"active,  t,  is  to  « 

teraperat;       tri  ■      *  sludin   b   r   Latio: 

is  resul      t   n  be  compared   it   the  zei  -  -      m 
solution  for  comparative  analysis  of  tne  effect  .on 

on  the  overall  proole  . 

ie  solution  for  this  n©v  torapcrature  distrJ         11 
require  a  relterat     procedure  in  order  to  ap  1,/  the 

i  ter    i  the  problem*   Is  a  first  est"    ,  the 
•  >n  bem  eri  ture  distribution  will  be  use;1. 
It  should  be  noted  that,  for  this  prol  ei  »  the  first 
half  3f  each  of  the  throe  te       ;c  distribution  equations 
}  and  (24)  will  not  be  affoctc-     )  ian       the 
•atur   istribu    .    :e^e  (q/A)  terras  may 
evaluated  once  and  for  all  at  the  cive  i  pro  >erties  a 
(q/A)  val   .      ictf  the  jvalua 

in  the  zero-radiation  a    ion,    lese  results  are  ta 
• 


The  3gc     art  of  each  of  these  equations,  involvi 
radiation  flux  (q/A)  ,  does  depend  on  the  temperature  dis- 
tribution,     terms  outside  the  integrals     be  evaluat 

at  the  :1vom  properties.  These  coof ficients  arc  a] 

tabulated  .1  '. 

Is  brings  up  the      roblom  of  numerical  evaluation 

of  radiation  (q/a)      ral  terms*  The  first  job  is  t 

r 

determination  of  (q/A)  as  a  funotion  of  radius  using  t 
selected  first  est      ;eiaperature  distribution. 

This  is  the  start  of  the  reiteration  procedure  to 
establish  a  final  compatible  temperature  distribution.  Each 
of  the  following  steps  will  be  repeated  until  that  temperature 

.3trlbutlon  is  reached,   ample  data  and  calculations  for 
each  of  these  3teps,  based  on  the  final  trial  te    \ture 
distribution,  are  she    \   Appendix  V. 


Steo  1 


4 

Calculate  the  '  «cfT  distribution  with  radiu3. 

Step  2 

Calculate  the  control  volume  total  radiation     sion 
flux  as  a  function  of  radius. 

(q/A)   ,   .   •  -  4(X)(KR)  (40) 

emission        S  6    1 

The  calculations  of  KJ'   ,    KJ    anc    (x)     are  shown  in 
Appendix  III. 


Step  3 

Calculate  the  wall  contribution  radiation  flux  ac  a 

function  of  radius  fro-  equation  (42).  These  results  will 

be  constant  for  this  particular  probls    ice  B  is  Hold 

w 

constant  as  a  given  boundary  condition* 

"teg  4 

Calculate  the  ;as  contribution  radiation  flux  as  a 
function  of  radius  from  equation  (41)*  Tain   will  require 
an  integration  for  each  radius  point  sinco  E_  is  in  t 
integral. 

Step  5 

Combi        to  determine  the  let  radiation  [*]    b  the 

control  volume: 

(o/A)„  ■  (q/A)  +  (q/A)   .-t   (q/A)   .   ,  (40) 

r       ga      ..^4.1      omission 


ot  the  results  of  stop  r-  on  a  large  3calo  plot  of 

Cq/A)  versus  radius  an.d  faro  in  a  curvo  to  determine  more 

complete  (q/A)  radial.  di3tribut3    [ata  -  ssnocially  near  the 
r 

wall, 

Step  7 

Calculate  the  (q/A)„  terms  inside  the  Integrals  of  the 
throe  temperature  distribution  equations  (20)#  (22)  and  (24). 


Step  8 


Plot  too  results  of  stop  7  versus  (y+)  and  integrate  these 
curves  in  steps  adequate  to  defino  the  temperature  distribution. 


r!ttl1        ■  suit  in  -  e      ate 

coefficients  previously  ei  "      '  and  tabulate         flix 

.  This  results  in  the  complete  "       o:1  t-  :  (  /  ) 

radiation  heat  flux  be:  erature  distribute 
equations  (20),  (22)   and  (34). 

Step  10 

Combj  e3 c  radiai  (o/A) 

tori:S  tabulated   in  Appendix 

itribution  with  radius  from  the  a.roro  >n: 

or   (   0  <y*  <  5)     use  equation   (20). 
For   \   5<y+<50)     use  equation   (22). 

For   (30  <  Sa*'  )     use  equation   (24). 

o 

ben  stops  arc  repeated  unt 
perature  distribution  la  nearly  Identical  with     tnter- 
temperature  distribution  of  that  trial •  This  particular 
sample  rro'  '  trial  s<     >ns  in  order  to 

roach  convergence    \  ill   Lus  cr  minus  one  decree  for  a 
tcept  th       kline  position.  The  final  temper 
Lstributlon  is  plotted  with  the  original  soro  trial  temp- 
erature distribution  for  comparison  In  1  ure  3C  . 

EHAS3  D  -  ITet  '-p. ■.Mat ion  Heat  ,;>ranrfcr  to  the  Y7aH 

The  next  step         a  trans:  sr  % 
calculation  of  the  overall  not  radiation  heat  flux    sb 
actually  reaches  tho  v;all  from  the  whole  cas.  !Tow  that  the 
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i]  oraturc  distribution  in  th    i  io  01    ished9  this  is 
a  relatively  si;     >rocoss. 

io  one  way  ,  as  rt    ion  anal       the  same  a.'> 

previously  evaluated  lor  the  wall  radiation  r     o  by  a 

£as  cube  in  steps  5  to  7  of  Phase  Bf  except  .for  some  changes 
in  subsorip     A   the  a    Ion  of  one  more  into;  ration  across 
us  of  the  ;  ■■*•  • 


Step  1 

valuate  tho  radiation  amission  (cTqL,  which  is  rocoived 
a  pipe  wall  area  element  cU  from  a  snail  cube  of  -;as  of 
dimensions  r.da  by  dr7  by  dx,  located  at  r-,  0,  and  x-,  • 


s. 


T(   r*+IT-  2r}   R   cor    ft  -f  x£    )   ' 
can  be  simplified  in  dlrionsionless  form  to  the  following: 


^  «   '    R.   R   '   R     PJtvX  ]fe 


Integrate  this  expression  along  the  complete  length  of 
the  pipe  to  obtain  tho  radiation  emission  (d''q)   which 
roaches  tho  wall  olomont  from  a  lor.;;  thin  rod  element  o 
located  at  T-y   and  &<*    •    e  axial  temperature  variation  is 


neglected  for  this  Integration.  Variables  are  chanced  to 
avoid  infinity* 


&* 


"/a 


(d°q) 


V 


ir' 


^        K     K     K  kJ 


vOo 


Ra  itf\e   "■?«*«*«,  (49) 


This  integral  has  '.    carried  out  as  previously  discussed. 

The  results  arc  shown  in  Figures  XI  throi 


Stop  3 


o 


atermlne  the  radiation     'ion  (d~q)  received 

Q 

by  the  T/all  element  fron  a  thin  cylindrical  gas  shell  located 

at  r-^  by  integrating  the  proceeding  e:        in  B. 

6  **$ 


(d2q) 


w, 


'/a 


\ 


^^/XM^W'-^^tele^^wyiW80* 


This  integration  lias  also  been  carried  out  as  discua    Ln 
step  7  of  Phase  2 •  The  results  are  shown  in  Figures  XVII  and 
XVIII. 


Step  4 

Integrate  this  last  expression  across  the  radius  in  order 

to  obtain  the  radiation  emission  (t3.q)r;  which  is  received  by 

tho  wall  element  fro::  the  whole  ,r;aa  in  the  pipe, 

r.  -- 1 


S   I' 


K 


4 


(a,l.^^^^^x^V'■5MW^*^cS)^,a, 


•  50  - 

The  radiation  emission  received  by  the  "wall  oloment  f*om  the 

gas  per  unit  of  wall  elomont  area  is  simply  (dq)w_  •  Since 

£ 

dA 

this  element  is  typical  of  all  such  v;all  elements  at  this  cross- 
section  then  (dq)VJ.  ■  (o/A)w  • 

This  result  should  also  be  corrected  by  the  real 

transmit tance  actor  (x)  •  This  loads  to  • 

6 

following  corrected  result. 


>  -.  ■   /i 


fc 


(q/A)w  = 

E 


^(wta*>  6-^-ok^^^4^^y^) 


7T 


This  Integral  has  not  been  shown  before  in  this  analysis. 

The  results  of  this  3    ration  are  sIior.ii  in  Appendix  VI • 

3 top  5 

The  second  component  contributing  (negatively)  bo  t] 
overall  heat  flux  reaching  the  wall  is  the  total  radiation 
emission  from  the  wall  element 

6 

Ginco  thG  "walls  are  assumed  to  be  "!  ,  €  s  i.  fJLso»  this 
wall  element  is  typical  of  all  such  elements  at  this  croc  re- 
section so 

(q/A)r  «  -  (54) 

6 


Stop  6 

How  conbiiio  terras  to  find  the  overall  not  radiation  heat 
flux  reaching  the  wall  from  the  whole  :-as, 

(q/A)   -   (q/A)w  +  (q/A)  (55) 

r        g      e 

The  numerical  results  of  this  analysis  are  shown  in 
Appendix  VI. 

IHASE  E  -  Total  ifet  T!oat  Transfer  Fate  To  The  '  all 

The  total  not  heat  transfer  to  the  wall  is  not  the  sasve 
as  the  net  total  heat  transfer  to  the  gas  control  volume 
located  at  the  wall*  This  is  due  to  the  difference  . 
radiation  influences  caused  by  the  angles  of  incidence  of 
radiation  beans  on  the  wall  area*   In  ad< it ion,  a     ;e  of 

a  occurs  for  the  convection  component  leaving  the  ^aa  and 
entering  the  wall*  This  causes  an  additive  effect  of  con- 
vection fro:::  the  local  -as  ring  control  volume  a      Latlon 
from  trie  whole  r.as* 

For  the  ;as  next  to  the  wall 

(q/A)Q  s  (q/A)0  4-  (q/A)p  (56) 

convection  radiation 

(q/A)     is  the  assumed  constant  heat  trans  7:  .be  for  the     . 

next  to  the  wall,    (q/A)0  s  -14,500  Btu/btr.ft 
(q/A),,  ,   at    (  ri    )   «  1,   lias  also  been  calculated  and  is 
tabulated   in  Appendix  V,    (q/A)r  s  +  2,158   BtuAr.ft*-   . 


Thorefore  (q/<\).  ,  the  convection  co-  >c 

c 

determined  from  equation  (56) 

(q/A)   ■   (q/A)   -  (q/A)  •  -  16,663  BtuA]  . 

C  O  2* 

Now  for  the  wall  heat  transfer   rate, 

(q/A)  «  (q/A)w  -f  (q/A)  (57) 


it  lii  this  case  the  convection  torn,  although  equal  in 
magnitude  to  tho  preceding  gas  convection  tern*  is  « 
opposite  Bign  (going  Into  the  wall). 


( q/  { )w  s  -  ( q/A )  c  =  -r  1  u ,  G  68  \ i  t  u/nr .  ft 
c 


.::  r.'all  radiatio  •   evaluated  in  Phase      • 

bs   arc  wn  in  Appendix   i    .      (q/A),,,  =—1667  Btu/for.ft 

Wp 

horefore»   the  total    .         .^at  transfer  rate  'jo  the  wall, 

incl-  T\     La t ion  for  tho  specified  problem  is 


(a/A)     ■     18,356  Btu/nr.ft 

m 
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10  combination  radiation  and  convection  h      xsfer 
problem  la  sc       considered,  for  approximation,  as  a  sun 
of  tvr>      ondent  solutions.   hen  the  temperature  distri- 
bution is  not  actually     rmined#    a  radiation  solution 
is  based  or:  a  constant  bulk  temperature  assi  "nod  t      -hole 
•as.  The  convection  solution  is  taken  from  the  results  of  a 
pure  convection,  zero  radiation  analysis  usin  ;  the  appropriate 
bulk  temperature  of  the  as. 

rJo  compare  the  sum  of  these  independent  a     Lmate 
solutions  with  the  combined  solution  numerical  result  s,  a 
oommon  bulk  temperature  of  the  :as  is  computed  by  eq  Latlon  (44) 

fxom  the  combined  solution  results.  T,  ■  1949°R, 

b 

Fluid  properties  are  assumed  to  have  the  same  value,  and 
all  boundary  conditions  are  the  sane.  Therefore,  the  pure 
convection  coefficients    heat  transfer  for  the  -as  to  the 
v;ali  Is  unchanged:  '.  '   10. G  Ltu^ir.ft  °P  Cro~:   equation  (<15). 

en  for  the  Independent  pure  oonvoction  heat  flux  t£  t] 
wall,   (q/\V,  ■  h('  b-T0)  ■+13,090  Itu/hr.ft2  . 

8   independent  constant  gas  temperature  a      ate 
radiation  contribi    .  to  the  wall     ~:>e   given  by: 
oAdams  6j  ,  paces  39  and  GO) 

(q/A)v/p   -  3  {  ^  -  x^T*  )   wheH6Tg  ■  Tb  ■  1949°R 

This  simple  calculation  yields  (q/A)r  -  -+  2127  J  tu/ir.ft 


-  .  •;  - 

Tho  sum  of  these  tv;o  independent  solutions  should  re- 
present an  approximate  total  not  heat  transfer  rate  to  tho 
wall, 

(q/A)  ■  (q/A)w  +  (q/A)„  s  10,117  Btu/hr.ft2 

For  comparison,  tho  oombined  analysis  equivalent  results 
are  given  as  follows:  {for  the  same  gas  bulk  temperature  and 
the  sane  wall  temperature). 


'qA)_T_+.,_  ■  <qA>    =  +  1^,668  Btu/nr.ft 


convection 


fo^MdUtlon      "  (q/A^H     ■+*•««    BtuAr.ft 


(q/i)  ■   (q/A)         ■  -+16,355     Btu/Vir.ft 

■  0  j     1  0  u^lJ.  »- 
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iSULTS 

and  T  isoussion  of  results 


results  of  this  thesis  may  be       J  as  follows: 

(A)  A  general  solution  for  the  -as  radiation  hoat 
transfer  problem  for  a  :ivon  gas  temperature 
distribution  in  a  eircular  plpo. 

(B)  Coneral  numerical  results  of  the  prol     ry 
integrations  r  iuiro<S  In  t       9  solu    ..br  a 
ran."o  of  problem  variables. 

(C)  A  procedure  for  Incorporate  Lon 
solution  with  the  convection  heat  transfer 
solution  in  "    to  ohtad      overall  combined 
heat  transfer  solution. 

(D)  a  s  3al  solution  for  the  combined 
radiation  and  convection  heat  transfer  proble 
by  the  methods  developed  above. 


(A)  The  eneral  solution  for  i;nc  gas  radiation  heat  transfer 
proble:..  is  .  iven  i     ail  in  Phase  ...  and  D  of  the  procedure. 
rho  only  restrictions  imposed  by  this  solution  are  as  follows: 
!♦  tsry  of  the  jas  anclosure  is  that  of  a  very 

io   oircular  pipe 
2.  The  radiation  properties  are  constant,  so  limiting 
the  range  of   F^as  temperature  extremes  in  any  .-iven 
•oblera 


P.   The  pipe  wall  Is  assumed  to  be  black. 

4»  The  temperature  distribution  is  specified. 

For  convenience,  but  not  of  necossitp,  the  solution  has 
been  written  with  the  additional  assumptions  that  the  tor.ro- 
erature  distribution  is  ar.isymetric,  and  that  the  tonperature 
ad  lent  in  the  axial  direction  is  relatively  snail. 

This  solution  provides  the  answer  for  t  .■'•.ion 

heat  transfer"  rate  to  any  selected  control  volume  element 
the  pas,  (as  shown  1  '..:\sc   '  of  the  Procedure)  includinG 
the  effects  of  all  b    is  in  the  pipe  and  the  entire  pipe 
wall  area.  Also,  the  net  radiation  heat  transfer  to  any 

lven  wall  area  element  from  all  the  par:  in  the  pipe  is  -ivon 
in  a  alr.iilar  expression  ('^lase  P  of  ."rocofiure }  • 

(B)  The  preliminary  integrations  required  in  the  above  sol- 
utions have  boon  carried  out  numerically  for  a  ran.  ;o  of 
problem  variables  generally  independent  of  any  s  A   led  tea  - 
eraturo  distribution.   In  partictilar,  the  first  integration 
alonp  the  length  of  the  pipe  has  "boon   carried  out  in  each  caso 

for  a  range  of  (KLR)  values  from  0  to  -s  10.  These  results  are 

a 

given  in  their  most  useful  form  in  the  oomi-lo  '  plots  against 
(H)  in  1 .Inures  VI,  VTIf  an  T.J  for  the  pas  control  volume 
cube,  and  in  Figures  XIII,  XIV,  XV  and  XVI  for  the  wall  area 
element.  These  results  may  bo  usod  over  a  vide  range  of  pipe 
radius  values  and  a  wide  range  of  gas  absorption  coefficients 
<Ka). 

The  second  integration  for  each  case  is  carried  out  in 
the  9  direction  around  the  pipe  at  constant  r.  r?herofore, 
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tho  tomperaturo  distribution  still  remains  outside  t" 

into~ral  for  tho  ax*    trie  caso.  These  into  "rations  aro 

carried  out  for  one  -    .0  value  of  the  dlnensionless  parameter 

(    ■  0.78).  The  results  arc     .  1  r  i  jures  X  an: 
a 

Tho  a     lea  ling  up  to  these  second        tio"is  are 
rather  invoi.    le   to  tho  roqai:        shiftily  t 
abscis aa  of  the  first  integration  result;:  bo  the       riate 

values.  Phase  B  of  re).   Also, 

the  added  complication  of  dealing  with  inf.  for 

several  of  the  curves  requires  tho  use  of  approximate  methods 
uutaerical  in  integrals  of  those  curves, 

other  integrations  in  each  case  across  tj         require 
tho  introduction  o±   tho  temperature  distribution  and  therefore 
c(m   not  be  yivon  in  advance  for  any  yenaral  case. 
(C*)  The  combined  convection  and  radiation  heat  transfer 
solution  is  dovelooed  in  general  fori,:  in  Phase  A  of  the 
Procedure.  This  solution  is  limited,  ai  shown  in  tho  develop- 

xt»  to  a  special  bu1    .rly  typical  caso.     particular,  it 
is  voted  hero  that  the  radiation  heat  transfer  at  the  wall  is 
assumed  to  be  relatively  small  compared  with  tho  convection 
component  of  the  overall  heat  transfer  at  the  wall. 

Tho  von  ICarman  const  ait  properties  pure  cc  it 

transfer  solution  10  v;as  used  as  a  basis  for  this  co  T.  iod 
analysis*  This  solution  wa3  chosen  as  tho  best  anal;  tic 
solution  for  the  convection  problem,  which  pormil    a  rel- 
atively simple  form  of  integration.  'Tic  Tar  tine  Hi  solution 
;   lich  covers  a  wider  range  0    >licationf  is  not  siy- 
"icantly  bettor  for  ^ases  (in  the  region  of  Jpj,^"  0.70), 
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but  requires  a  m  ore  c  of  in'     tie  .   lie 

temperature  distribution  givon  by  the  methods  of  the 
voi  Karman  solution  is  not  as  reliable  near  the  contorline 
of  the  pi  e,  but  thli    akness  was  not  consider       Vficant 
enough  to  require  the  use  of  the  'lartinclli  soluti    s  a 
basis  for  the  combined  analysis. 

The  blgges  the  combined  analysis  concerns  the 

evaluation  of  the  (q/A)  radiation  heat  flux  term  as 
component  of  the  tota:\  heat  transfer  rate  at  any  poinl  .    ;re 
appear  to  bo  two  possible  interpretations  as  to  tlio  mea 
of  the  apparent  heat  flux  at  any  .riven  point.  Those  two  con- 
cepts are  as  follows: 

(1)  The  apparent  heat  flux  (q./A)    rep         ic  total 
heat  transfer  rate  c  constant  radius  boundary 

of  a  given  rj  . -■  control  volume.    tat  is,  it  represents 
the  convection  heat  transfer  rat    r  unit  area  at  that 
boun«"    [at  r  ),  plus  tb    >tal  let  contr:   ion  of  all 
the  radiatic     it  transfer  rate  from  all  inner  control 
volumes  (at  r  :  r-j)  which  actually  crosses  that  boundary, 
per  unit  boundary  area,  as  radiation  heat  tra     ■   his 
ic  the  most  strai;  htforword  anr1  phyal 

interpretation,  but  unfortunately  this  r        r  lation 
function, as  s  oc"    '-,1s  very  difficult  to  obtain, 
di'    ilty  Is  brou  vvt  about  because  of  the  required 
account in-         transmittance  of  the  pas  and  for  the 
adjustment  of  the  contributions  of  radiation  from  the  gas 
seated  upstrea       nstream. 
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(2)  The  second  interpretation  defines  the  apparent  heat 
trans    'lure  (q/A)    at  any  point  to  bo  a  net  heat 

clpp 

transfer  rate  ^or  the  'as  located  at  that  point*  That 
is j  it  represents  the  algebraic  sua  of  the  eonveetic 
hoat  transfer  rate  per  unit  area  at  that  point,  {at  r_), 
plus  the  net  radiation  hoat  transfer  rate  per  unit  area 
(normal  to  the  wall)  at  that  same  point.  This  statement 
Is  not  complete j  however ,  until  some  depth  is  assigned 
to  the  ;/;as  under  consideration  at  this      ,  since  the 

as  radiation  must  refer  to  a  volune.    Ls  is  accomplish- 
cd  by  referring  the  radiation  heat  transfer.  Cq/A)p,  to 

a    ■' 


*-v 


tensionless  unit  thickness  of  constant  t     ature 


as  located  at  tills  point  (at  r,  ).  ?hia  is  the  (q/A) 
tern  which  is  used  in  this  development ,  as  defined  at 
the  end  of  Riaso  P   of  the  .Procedure. 


It  should  bo  noted  at  this  point  that  t"         nces  in 

interpretation  of  the  (<l/A)-__  and  (q/A)  terms  have  no  effect 

app         i» 

on  the  developed  temperature  distribution  equations  (20),  (22) 

and  (24),  except  insofar  as  ishat  is  now  raoant  by  (q/ )     = 

(q/A)    for  the  j;as  at  the  wall.   If  (q/A)  _  is  ta]     i  the 
app  app 

first  uense,  then  (q/A)n  ■  -  (q/A)  .  That  isj  the  total  c  crib  in- 

eel  heat  flux  for  the  ;as  at  the  wall  is  equal  in  magnitude  to 

the  total  combined  heat  flux  to  the  wall  (as  In  the  pure 

convection  case).  Dut  if  (q/A)    is  taken  in  the  second 

app 

sense,  then  (q/A)Q  f   -(q/A)v;.  That  is,  the  not  combined  hoat 
flux  for  the  gas  at  the  wall  is  not  of  the  same  magnitude  as 
the  net  oorabined  heat  flux  to  the  wall.  This  second  case  is 
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the  case  use         \  this  development*     The  re]  ishlp  between 

(q/A)     and    (q./a)     Tor  tills  oaso   is    ;iven  in  of  the 

o  w 

Procedure. 

("' )  In  order  to  demonstrate  the  application  of  the  combined 
solution,  a  numerical  sanple  problem  was  carrier!  to  completion 
for  pure  carbon  diozlde  gas  a3  discussed  in  Phase  C,  Phase 
and  Phase  E  of  tho  Procedure. 

As  a  starting  point,  a  pure  convection  solution  was  car- 
ried out  to  determine  a   '  Ltial  trial  temperature  distribution 
(Figure  XIX),  and  to  establish  boundary  eondil   a  for  the 
combined  convection  and  radiation  problem. 

The  combined  solution  the      "red  a  re    'ation  pro- 
cedure  to  establish  a  compatible  temperature  distribution, 
is  process  is  described  in  ?haso  C  or  the  Proceduref  and 
sample  calculations  are  given  in  Appendix  V.      resulting 
temperature  distribution  is  :;Iven  in  "l  ~utg   -.  . 

.eneral  o'liape  of  the  ocmbined  solution  temperature 
stribution  curve  is  si  li  ioantly  different  from  the  original 
zero  radiation  result •        aperature  appears  to  come  closer 
to  a  uniform  value  >e,  so  that  the  temperature 

gradient  near  tlae  wall  is  relative-.    oator  Tor  a  given 

turn  temperature  or  bulk  temperature  of  bhe  as. 

ls  drop  off  oi'  temperature  toward  tho  center  of  the  Pipe 
is  of  special  interest.  This  is  cue  to  the  relatively  favor- 
able position  of  this  us  for  radiation  interchange  with  the 
rest  of  the  -as  in  bhe  1  c,  coupled  with  the  earlj-  assumption 
that  (o/A)a"  varies  linearly  with  bhe  radius.    3  integration 
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procedure  used  In  the  v    iraan  method  of  analysis,  as 

extended   to  this  application,  provide  a  relatively  weak 

•oorature  distribution  result  near  the  eenterline  (as 

d  by  the  .       "iation  result),   Tb'.;ov^,  it  la 
believed  that  the  drop  off  In  te      re  near  tho  center- 

s  in  fact         en  the  a::ial  temper  itva?  '   »nt 

la  nearl  !  is . 

he  :.     'ature  i    Lculation  rj       the  con 
not  yield  a  o    '  IbV-     irature.   his    nperati        at 
si  :ilf  icaiitly  Influence  the  sther  temperature  calculations  be- 
sero  area  asse elated    h   the  centerline  (r  s  0). 
Therefore,  It  is  c  .  ' "  -  »d  justifiable  te  faro  in  an 
approximate  temperature  for  this  singular  point.    icr  nore 

methods  of  intc  -ration,  such  a  lolli 

e  convection  analysis,  might  t       to  obtain  -     accurate 
*y  near  the  eenterline.  However,  the  influence 
or  t    a      MP  '         near  the  ce:  1      ■  is  very  small 
in  the  >verall  heat  transfer  a     is, 

t  bal  net  heat  tra     ■  rate  to  t        is  oal- 
oulatc   "   i.as^s  ]  and    )f  the  Pr;  celurc.     i  result 

•  "  bhat  bh  radiation  effects  led  to  a  ;  bal 

it  transfer  rate  to  the  wall  Tor  the  iv        ary 
mst  The  tern  >erature  distribution  curve,         ., 
temperature  of  the  gas  is  lower  foi 
i>ined  solutJ         >•  Therefore,  tho  conclusion 

he  overall     st  of  tho  rad     i,  3 
bh  a  zero-::-  ,  is  to  increase  the  rate  of  total 


not  heat  traj        t       for  a      bulk        bure  o 
hot  gas  *     ild  . 

:  concl  i;      '   must  be     lasiso  ,  ^n  a 

atively  restricted  solution  reach  the 

Ltatlons  and  assumptions  previously  o       .       be- 
lie   that  the  :  ioa '  s<  '   3  If      -  -   Imposed 

ted  heat  t]    or  eoluti   I   bhi  re  triction  to  constant 
fin    roporties  thron  "         a.]  sis.  >f 

riable  propertios  may,  in  Peat,   be  very  si  ;  "      ; 
especially  in  the  bound    layer  ire 

nt  is        .-  >.   It  is  still  wssible         o  "-".nod 
°ects  of  r  13  at  '  •   "        "  o 

properties  co  .    red,  mi  nt  ov  blned 

analysis  to  a  reductjl   in  the  total      eat  ti 
to  the  wall  in  a    oases  '        "   this 

constant  proportios 

'ore  conclusive  r     is  arc  is 

combined  analytic  so]     i  with  the  su  ndent 

conveotion  and  simple  x    ition  solutj  re- 

qulre  complete  to         diotributi  .  To 

accomplish  this  comparison,  e  bulk  temperature  of 

the  yao  v.ras  hole    stant  for  the  two  cases  as  described  in 
Phase  ';  of  the  Procedure.  The:     total  :  at  t! 

wall  was  compared  for  the  two  methods  of  solution* 

This  co  tparison  s    "  that  the  total  net  heat  transfer 
rate  to  the  wall  was  higher  In  the  combined  solution  than  in 
the  summation  of  bl    dependent  solutions*        I  lion 
component  of  the  heat  transfer  was  lower  at  the  wall  in  tl 


-•       OxJ 


CO 

transfea  ,      >  sfor 

rate   to    :  1  '  Pc  :  tlio  cc  6   ana  ,  >lutlon« 

arefor    , 

solui  *0      "^10 

*  latii  *    felon  e  tion  heat  .cnu 

jod   in  a  c         i         a  a]  jrature 

ti  ■      '  •'■        "  sly  de   aril 

transfer  soluti  -  l« 
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qoTJCLUsio 

Then;ial  radiation  snorts  a  significant  influence  on 
the  total  net  hc.it  transfer  <.':  ■     turbulently  flowing  hot 
fluid  In  a  circular  pipe*  A  s       to  this  c 
radiation  and  convection  heat  transfer  ^robloK.  is  presented* 
Constant  flui     >perties  arc  assi     throu  hout  the 
analysis,  A  sa:.i;jlo  numerical  probl   is  carried  to  com- 
;  for  pure  carbon  dioxide  ..-as. 

The  results  of  this  constant  properties  solution 
reveal  e      Meant  gain  In  the  total  net  heat  transfer 
to  the  wall  eonpared.  to  an  equivalent  sero-i'adiati,     lution, 

le  simple  addition  of  lndope     i  radiat:     i" 
convection  solutions  1  ilk  temper      do  3 

not  provide  an  adequate  total  net  heat  transfer  solution, 
!'hc   couplin  ;  effects  of  radiation  and  convection  oust  be 
introduced  in  a  corl>inod  analysis  to  determine  the  tornp- 
eraturo  distribution  in  order     'equate!    ,cribo  the 
heat  transfer  solution. 
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Tho  results  of  this  thesis  cannot  be  considered  an  a 
conclusive  explanation  for  the  effects  of  rad'atiou  on  the 
general  problem  of  heat  transfer  from  a  turbulently  flow! 
fluid  in  a  lciy  circular*  pipe,   Accord;   ]  ,  tho  following 
recommendations  are  presented  as  the  most  promisin  ;  ave 
for  significant  improve    :  in  the  understanding  of  this 
subject. 

A.  Explore  the  feasibility  of  usl  -  the  Martlnelli  [5j 
forr.  of  analysis  for  the  temperature  distribution  equations 
for  the  co-vibiiied  radiation  and  convection  analysis. 

added  refinements  of  a  f!artinelli  type  solution  woul   -rovide 

fcter  temperature  distribution  Information  near  the  center 
of  the  pipe.  The  effect  of  this  snail  center  section 

irature  distribution  on  the  ovorall  analysis  i3  relatively 
11  in  ar  case,  but  it  is   certainly  more  significant 
radiation  is  present,  "horefore,  it  is  believed  that  such  an 
iprovement  should  be  3    bigated.  The  significance  of  this 
>rover".ent  would  be  raised  for  cases  in'  whioh  radiation 
represents  a  higher  contribution  to  the  overall      transfer 

B.  Extend  the  combination  radiation  and  convection  heat 
transfer  analysis  to  include  the         f  variable  properties 
of  the  fluid.  It  is         that  a  si  nif leant  Improvement 
would  be  obtained  by  intr       variable  properties  for  the 
convection  properties  of  the  fluid.  To  include  variable 


radiation  properties  of  the  fluid  would  require  a  much  more 
co:     radiation  solution,  but  the  improvement  here  wou    rob- 
ably  not  be  nearly  so  significant  for  most  problems. 

Therefore,  it  is  rec    ruled  that  only  the  properties 
associated  with  the  convection  heat  transfer  b     sidored 
riable.  'This  will  lead  to  a  velocity  distribute    hieh  is 
unction  of  the  temperature  distribution.    b  effect  of 
this  consideration  will  be  especially  apparent  in  the  boundary 
layer  region.  The  overall  solution  may  yield  significantly 
different  results  -  especially  for  larrer  temperature  extremes. 

C.  Carry  out  the  combined  n  diation  and  convscti 
analysis  using  the  total  heat  flux  concept  for  (<j/a)ft  _»  the 

apparent  heat  flux  at  any  point,  instead  o  at  flux 

concept  used  in  this  analysis.      definitions  of  these  two 
concepts  arc  disci  bhe   ssults*     ause 

of  the  more  tangible  physical  picture  offered  by  the  total 

it  flux  concept ,  such  a  result  could  be  rr.orc   oasiT'y  inter- 
preted,  and  would  therefore     ide  a  stronger  a      b  for 
the  solution.    e  objective  of  this  new  more  c       .nalysis 

",.ie  to  provide  another  numerical  solution  for  a  check 
with  the  results  obtained  usinr  the  analytically  simpler  net 
heat  flux  concept. 

The  chief  probl    iich  is  brought  ux>   by  this  new  anal;  - 
sis  is  the  evaluation  of  the  radiation  component  of  tl 
apparent  heat  flux  at  a  .;   ivc       •  At  first  thought ,  one 
might  suspect  that  the  total  net  radiation  flux  croc  '    a 

■/en  boundary  at  r.,,  is  simply  the  sum  of  all  net  radiation 
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effluxes  from  all  the  ring  control  volumes  located  at  r  (  r.  , 
summed  from  r  »  0  to  r  •  r,  •  There  arc-  two  things  wrong  with 
this  simple  summation.  Firstf  and  most  obvious,  axial  radia- 
tion components  are  not  directly  accounted  for  in  this 
summation.  Fortunately,  this  difficulty  is  not  s ifjnif ioant 
for  the  infinitely  long  axisy metric  pipe,  if  the  axial  to;  - 
erature  gradient  is  assumed  to  be  very  small.   In  that  caoo, 
mutual  syraetry  exists  for  all  control  volumes  alon£  the  axis 
at  a  given  radius;  such  that  all  the  net  radiation  efflux 
from  a  series  of  control  volumes  at  a  -iven  cross-section  may 
be  assumed  to  flov;  in  a  radial  direction.  Therefore^  the 
axial  c:    tents  arc  accounted  for  in  this  indirect  manner. 

But  the  second  problem  is  not  conquered  so  easily.  This 
difficulty  lies  in  the  coupling  effect  of  radiation  a 
convection.  The  sum  of  the  total  not  radiation  heat  effluxes 
from  all  the  control  volumes  inside  of  the  radius  r.,  does  not 
necessarily  cross  the  radius  r^  as  radiation  heat  transfer. 
The  mechanism  of  heat  transfer  can  change  to  convection  heat 
transfer  for  a  part  of  that  net  radiation  efflux  wi  ich  is 
absorbed  by  the  interim  ~as  before  reaching  the  boundary  at 
Tn*     Therefore,  the  sum  of  all  the  radiation  heat  effluxes  from 
all  the  control  volumes  inside  of  the  radius  r-,  represents 
the  total  net  radiation  flux  at  th.  boundary  at  r^,  plus  a  part 

of  the  convection  heat  flux  at  the  boundary  at  r  .   If  this 

1 

sum  is  then  added  to  the  local  radial  convection  heat  transfer 
established  by  the  temperature  gradient  at  r., ,  this  total  will 
represent  moro  than  the  actual  total  heat  flux  since  a  part 
of  the  convection  heat  transfer  is  accounted  for  tnice 
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A  simple  solution  to  this  difficulty  would  bo  to  1  ;nore 

it.   v-on  tho  now  radiation  term  (q/A)  at  the  boundary  r 

™  1 

would  bo  simply  tho  following  into,  ral. 

— 
(q/A)   -   r  (dq  )   dr 

O 

This  would  certainly  lead  to  an  approximate  rosu:t,  but 

its  sirriif  icaoce  would  bo  doubtful  until  it  is  compar    Lth 
re  enact  solutions. 

The  more  exact  solution  for  this  expressira  could  be 
solved  by  the  folio;..    jroceduro: 

Carry   out  a  complete  three  dimensional  analysis  for 
each  control  volume  cube  at  r  inside  tho  boundary  r., ,  but 
only  compute  the  net  radiation  which  actually  penetrates  the 
boundary  radius  r-,.     b  is,  take  into  account  all   .  -way 
wall  radiubion  which  reaches  the  cube,  plus  all  one—way  ^as 
radiation  from  the  pas  looat     stween  rn  an;"  r  which  reaches 

JL  O 

the  cube,  plus  the  fraction  of  the  total  omission  from  tiie 
cube  which  actually  crosses  the  radius  r.  •  This  last  emission 
term  amy  require  a  rather  involved  calculation  to  account  for 

trans mittanco  through  sono  average  path  length  to  the 
boundary  at  r,.  Tho  net  result  of  these  integrations  would 
yield  a  modified  heat  flux  tern  for  that  cube  (  da    )r   for 
a  specified  value  of  the  boundary  radius  r  •  'Then,   this  sane 
evaluation  should  be  repeated  for  other  representative  cubes 
inside  r-,,   bach  cube  will  bo  representative  of  ail  cubes  at 
the  sane  radius  at  that  cross -sect ion,  so  r     eat in  -  a  r 
control  voluiae,    -:t,  su     all  control  .        tributiona 
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to  tho  flux  at  tho  boundary-  r_  • 

Then  for  this  analysis: 

rl 
(q/A)  ■   r  {  da  )   £r 

o       -1 

?/hore  the  tern3  Inside  tho  integral  were  evaluated  for  their 
contributions  at  r  ♦  Therefore,  this  whole  process  mu3t  be 
repeated  for  each  selected  value  of  tho  boundary  radius  r  • 
This  would  finally  lead  to  a  (q/A)^  distribution  with 
radius  which  could  be  used  in  the  total  heat  transfer  rate 
concert  definition  of  (o/a) 

Then  a  reiteration  procedure  would  be  carried  out  to 
establish  the  compatible  temperature  distribution,  '.'his 
reiteration  procedure  would  require  recalculation  of  each 


?P* 


(<l/A)r  value  for  each  trial  because  of  the  variable  tenp- 
orsture  gas  contributions.  This  orocess  should  th 
converge  to  the  final  temperature  distribution  and  to^al 
net  heat  transfer  solution. 

Needless  to  say,  the  work  involved  for  this  job  would 
be  increased  about  tenfold  over  the  work  required  for  the 
solution  presented  in  this  thesis  (even  though  constant 
properties  are  maintained)*  However •  the  radiation  solution 
in  this  thosis  could  be  used  to  advantage  to  reduce  this 
work  load,  especially  for  the  axial  and  ->   direction 
integrations* 

It  would  be  interesting  to  seo  if  this  solution  would 
yield  significant!,    re  accurate  results  than  the  other  more 
simple  analytic  solutions,   \t  least  one  or  more  sets  of 


corporis  on  solutions  by  different  methods  will  bo  required 

In  order  to  show  which  method  is  most  suitable  to  provide 

significant  results  with  a  practical  anoint  of  work* 

A  final  extension  of  this  radiation  tern,  (q/A)„  , 

rl 

evaluation  could  be  ma    !ter  the  accurate  solutions  wore 
determined  for  several  cases.   Frobably  a  relatively  simple 
approximate  two-dimensional  breakdown  could  bo  developed  based 
on  the  results  of  the  more  accurate  three-dimensional 
solutions,  Such  a  result  v;ould  bo  valuable  for  general  engi- 
neering application  for  combination  radiation  an  convection 
heat  transfer  problems  in  Ion"  circular  pipes. 

A  final  rec emendation  concerns  the  neod  for 
experimental  evidence.  Accurate  experimental  evidence  is 
needed  concerning  the  effect  of  radiation  on  the  overall  not 
heat  transfer  rate  from  a  hot  ,.~a3  flowing  in  a  cool  pipe.   In 
particular 9  temperature  distribution  measurements  would  be 
moot  valuable.  In  this  wayg  a  strong  check  point  could  be 
made  to  compare  the  analytic  development  with  experimental 
evidence. 

Also,  comparative  experimental  evidence  should  be  obtained 
to  show  just  how  significant  the  radiation  effects  become  with- 
in the  range  of  practical  application  of  fluid  flow  heat 
transfer  proble     lere  convection  usually  dominates  the 
heat  transfer  process.  Or  from  the  other  extreme,  how  is  a 
radiation  dominated  heat  transfer  problem  chanced  by  the 
coupling  effects  of  added  convection  heat  transfer? 
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A  fev;  careful  experiments  -VI  n  ;  these  lines  should 
L dance  to  ,      tore  It  is  most  profitable  to 
aako  uso  of  combination  convection  and  radiation  heat 
transfer*  This  shou]  3  pr  vide  aii  answer  to   the  choice 
iting  :/^n  and  a  .       aon-radiatini 

pplicj  -ion. 


A   P   P 
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Idor  control  vol  ;   as  shown  1.       "  uros        ~. 
-  flor;  throu  ;h  this  control  volume  is  asc 
.  o'ic  ti        flow,  .ill  fluid 

constant*      ~ral  references :  Si  ,   •■   $        ,        , 
• ) 

1.  Contlnult  on 

; :  r  -'  r 

2.  'bstentun         ation  -  x  direction 


r>ee  Pi  otto 

a. 

(negleotin  -  all   second  order  derivatives    'n  the 


-^^H^W)  -fs^  ri  "*%?     (A-2) 


axial  direction  an  be  in"  relativ  11) 

Z*     Combine  equatJ  '    -1)  an      '.  -2)  by  a 


4«  .    last  tern  of  :i  ret 


•aclient  of  tin  lux  ii  action. 

n::all  c  the  other  terms 

;leeted. 
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5.  Introduce  the  timo  average  velocity  concept  for 
turbulent  flo:  . 


r 


V   «  V  -f  V 
X     XX 


v  s  v  -t-v 
r    r   r 


3-fV 


*v 


G.  Enter  these  terns  In  equation  (A-4)     retain  only 
oho  finite  time  average  terras  t  i  )  tain  the  following: 


I 
-2- 


dx    t  or 


Ut^X   •  jO       V1  v.'.  ) 


r  x 


d: 


6, 


(A-G) 


B.         ...... 

Consldor  the  control  volume  ao   shown  in  Figures   I  and  II. 
(Genoral  rof  »s:      2,4,7,,     and     0    ) 

1.     Apply   L.ie  first  lav/  of  thermodynamics 


-  dW       *     di  4-  d(    v^    ) 


C-l) 


ZfAr&rpA'f 


:    -2) 


whero  (  dg  )  is  the  not  radiation  heat  transfer  rate  to  the 


control  volu-  ait  volurjc 


Y         r  ^ 


r" 


rr\T 


(1-3) 


jl  +  <«  K>\1  :[r^+  ^  l^r^+^p^  ' 
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ote  the  continuity  equation   (A-l)    In  equation 


(.,-4)   and   simplify  to  obtain 


&.+&Q- 


<  \ 

o1  r     j 


^ 
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£.     Combine   aquations    ('■'-£),    (B-3)   and    (  -43  .  iclc 

by  2'0rdrclxdt. 


*    -• 


r^~%'  r^^-t +pn- 


N. 


^rf 


2.    Ac-   /On, 


if 

T 


^v-    A 
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4.  Combine  with  t  equation  ('.-":} 

simplify  to  obtain  the  following: 


\J-    (y±L 


d?r 


5.  For  pur  poo  oa  of  this  heat  transfer  anal-sis,  the 
pressure  gradient  term  and  the  viscous  work  borm  are  small 
compared  to  the  other  t     in  the  quation,  so  the;;  will 


be  no  ;lectccU 
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6.   tow  combine  with  the  continuity  equation  (A-l), 
to  obtain  the  followingi 


(  -10) 
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7.  toratu  j  axial  diroct- 

ion  is  relatively  small,  so  no /loot  t 

8«  :  or  introduce  the  tl         concept  for  turbulent 
flou  into  equation  (B-10)  velocity  and  t    ?ature  tori.13. 
tain  all  finite  tl  is. 

It  should  be  noted  that  in  waiting  to  introduce  tho 
tiro  average  conoe    \   tills  poI  it,  one  relatively  import- 
ant torn  la  not  3-,  ■  is  an 

-or  a  J  -'  xj  vV)  ^T::  /     ,  which  ml  ;ht 

I   t;     r  *  lr  J 

bo  significant  in  some  analyses.  Fortunately ,  this  ter:.i  is 
not  significant  compared  to  the  othor  remalninc  terras  in  tho 
energy  equation  for  purposes  of  this-  heat  transfer  analysis, 
and  its  omission  is  Justified.   his  is  due  to  t      h 
temperature  gradients  and  temperature  level  encountered  in 
this  application. 
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-.loral  references:  [1]  ,  [  [0]  .) 

The  oiiilsaion  coefficient    (11  )  and  the  absorption 

e 

coefficient  (K  )  aro  evaluated  as  functions  of  the  ;c.n 
a 

Lty  ancl  {*\S-)  >s«  The  dimension  of  each  or  t':C3o 

coefficients  is  the  reciprocal  of  length.  It  is  assumed  in 
thia  analysis  that  these  coefficients  are  equal  constants  to 
bo  evaluatod  at  sone  given  bulk  ten  >erature  of  tho  gas* 

KG  -  .c  *  K  (CI) 

€ ,.  s  (rO,.  (  1  -  o  "KL  )  (C-2) 


where  L  is  the  moan        th 

€^  ia  the  ras  c      .ty  at  the  an  bulk  temperature 
;ra;-c  (PL)  value 
and  (x)^  is  the  real  :a3  traus.iittancG  weijhtin  ;  factor 
described  in  stop  9  on 


6     Kg 


Solving  for  It 


=  -  X  /n  (1- | 


I       "  U)  (C-4) 


is  can  bo  simplified  to 


K   ■  i  »  (    6n   )  (C-5) 


'u   e 


-    7      - 


Thsn  solvirv;  for  tho  di  Less  .1  •.-', 


KR      * 


P  &  ( 

Z 


(C-G) 


B.     GEOMETRY  OP  T2I  IT  -  G. 

8n  Figure  XXI 

the  following  dimensions 
Angle  CDS    «  e0 
;    ■  r, 
AB     »  27>  ■  r 
AD     =  DC   ■  xQ 


CE     a  , 


'/     R^t  r?  •  2  Rroos0o 


AS     -,       r2.X2 


TV"? 


V 


'  R2-+  r2  -  2  PrcosP   +  x2 
1  o       o 


Ancle  ABE     =  3 

cos/3  =  cos   (A    .)  ■  (& )2-f-(r :\)g-  (A.:;)g 

2  (AJ  ) (I    ;) 


(-1) 


r>       <*i       o 


oos/9«  3- %.■■%. i^  -  a  V-*o-  ri 


2R  v'-    4  *?-  2Rrcos9  +  x 


(D-2) 


eos/?S        *"  rl  corj 


(D-S) 
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Figure  XXI 


cue*  fij 


A&^z  **&/?  ^t^s+^v 
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APPirmix  in 

TABULAT3  ■  ICAL  RliSUI/ 

A.   Convection  Analysis 

All  properties  oro  ovalucte^  at  the  gas  film  temperature, 

(reference  McAdama  [6] )  for     CARBON  DIOXIDE  OAS 


PARAMETER 

sy.: 

valui 

u  r 

Film  temperature 

\ 

1275 

OR 

Gas  pressure 

r 

1 

io  sphere 

Thermal  conductivity 

k 

0.0504 

BtuAr.ft2  °P/ft 

Specific  heat 

°P 

0.2 

Dtu/lb  °" 

Viscosity 

M- 

0.0311 

CenLipoises 

Prandtl  'lumber  «  c^V^ 

"ft 

0.70 

Reynolds  Number  »a' 

e 

20,000 

v* 

Pipe  diainetor 

D 

.1607 

foot 

Gas  density 

/* 

0.0473 

lb/ft 

Mean  gas  voloci 

V 

cs.o 

ft /sec. 

Friction  factor 

f 

0.0065 



vail   shear  stress 

T 

0.0134 

o 

/ft" 

Kinematic  viscosity*  V 

7 

1.591 

'ir 

fi 

Stanton  "umber  »       h 

"st 

D. 00391 

..... 

V* 

Hues  el  t  Number  •  hP 

IT 

54. S 

C oof fie lent  of  heat  t 

h 

10.0 

ntu/hr.ft2  °P 

fer   (zoro  radiation} 

Total  net  heat  transfer 

(q/A) 

-14,500 

3  tu/hr.ft2 

flux  at  wall 


- 


isionless  coordinate 

ratio    (;    /  =  1  ;  ojo      ) 

P 

a  / 

(assumed) 

C  o  iv  s  o  t  Ion  t  ornoratu" 

^(q/A)^ (V 


VA 


»  „ 


o 


1 


-    . 


1/ 


o 


;1  on   fcnal   s 

'  11  pr~  as  ter.ip- 

-t.  ,    (Fef. 

Braturc 

TQ 

ter 

Avon  ='  ••- 

*blo 

at   (r  x) 
at  2( 

radius 


Q.e3s  oalsa 
coefficient 


Lnensionleas  a    .       >tion 

co officio 


transmittonce 
v;c '        actor 

stent 


)  for 

- 

• 

yT 

' 

X5 

"! 

.: 

o.: 

L) 

%": 

ft- 

:  ) 

vj#  0  vjs^ 

c 

o.c 

€«- 

0#107 

— -— 

0,0853 

foot 

II  ] 

c 

0.78 

D.78 



(x) 


0.1130 


0,1713  x     va/ir.rtJ(°:  ) 


10' 


■ 


IV 


.,ori     adiation.    Pur  KJtlon  Te  Lstributlon 

Calculation 


Zona  I,    (  0  ^yV  5),   equation  {20},     I"-,  (q/AK  ■  -    J 


1 


(y/po» 


o.cx 


0 

1 

2 

!■- 

4 
5 


(otf^,)/      (TQ-V)«f1(q/A)0::  Col.   S 

0 


0.70 
1. 

2.10 

• 
3  •  '■-.' 


0 

-70 
-130 

-549 


G20 
509 

759 

820 


2one  IIJ    (5<y*-<30),   equation  (22),     *V(q/A)     •  -"S5 


(y/r0) 


y*"  5/n(l-x.*7;^KIY  jt)      !?,.-T)San  aboV( 


A 


0.0 


).05 


vJ 


10 

12 
1G 

20 

Of 


28.45 


0 

0.6 

X.  ' 

'. 

.116 
4.661 

. 
6.675 
7.2G5 
7.E 


0 

-35 

-175 

-.340 
-464 
-565 
-GG5 
-724 
-749 


1074 

1239 

1523 


5ono  1X7,    (3 0<yt<  r J ~~TT  equation   (24)     fTTqJ 


-   •.  • 


0.0527 

.15 
0.2 
0.25 

0.; 

0.35 
0.40 
0.45 
0.50 
0.60 
.70 
0.G0 
0.90 
1.00 


GO 
56. 

0i.     • 
11.     . 

142.1 

17C. 

109.0 

220 

2S6 

204 

542 

398 

455 

512 

569 


2.rA     +/50) 


0 

1.605 
.610 

4.729 
.07. 

.320 

. 

c.- 

.798 

7.557 


0 

-:. 

-260 
-3 


-471 

- 

— oo4 

-560 
-GOG 

- 

-706 
-755 


s  as  above 


1908 


2110 
2153 

21G2 
2203 

2525 
2554 


DIstrib  jefficien 

The  following  temperature  n  coefficients  for 

i  Int  •  ),    (22)   and 

based  ou  film  temperature  .  babulated 

Append'.:: 


diation  te^oraturc  VAJ/'*  TTS 

distribution  coeTTTc  b  : 


Zone   T,(0<-/<5)  ^Pr  0.00481     °P/  '.ft' 


/?c  oi/gor9p 

^     PV      /> 


Zone   II,    (r<    -<  30)  S  0.(  2     °  /  >tu/hr.ft2 

,  /  e  6o 


/' 


Zona   III,    (30<yr<   r*")  2. 0.017!         °p/  bu/hr.ftS 

/-  C..< 


/? 
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SA 


Step  1     -     Calculation  of 


rxA 


Assuricc" 


0 
0.10 

0.20 
C.30 
0.40 
0.50 
0.G0 
0.70 
0.80 
0.90 
0.05 
1.00 


(T/lOO) 


1950 

144,590 

1991 

157,141 

2010 

1GG,170 

2050 

169 , 818 

2030 

169, 

2020 

166, 

2001 

150,320 

1967 

149,699 

1012 

133,645 

179G 

104,510 

1651 

74,300 

550 

915 

,4 


24,7G8 
0  ,918 
2  ,465 
2?, 090 
29,090 

a  » 

27,463 
25,643 


17,903 

N   720 
157 


1r>     ryrjr-i 


W 


>tcp 


-     Calculation  of  Radiation  yjriission  fro.-?:  a  ' '"as  Cube 


Located  at   (r_F  )   , 


(qA),  ■  -4(x)  ay  )    , 


(x)      S  0.1138      ,    (K      )   ■  0.78 

c 


\/P 


E 


•(q/A) 


0 
0.10 

0.30 
0.40 
0.50 
0.G0 
0.70 
0.90 
0.90 
0.95 
1.00 


24,768 
26,918 
2  ,465 
B9f090 
29,090 
28,521 
27,463 

or;    oa  <z 

<_<*_• ,  ij  »«^ 

22,393 

17,903 

10,723 

157 


3794 
9557 
101 
10520 
10329 
10127 

9-; 

0105 


G357 

4519 

56 


- 


Ste.^  3     -     Calculatio  o  '  ay     a   Lation  rrow  '."holo 

ill  to  a     aa  Cube  r.-ocateO  at    (r  /■■)• 

1 

["Frora  Equation  (42)j»    (See  Plr^urc 


rxA? 


0 
0.1 

0.2 

>.7 

. 
0.0 
0.95 

:.. 


irr* 


(I  -  tl  <:»<. 


fel 


(%)^^M\(GU,Z) 


1.1330 
1 .  X 

1.1430 
1.1675 

i      .i.CL 

1.2733 

j.  .  < . 

1.4214 

1.5125 

1.7 

1.1400 

0.5313 


20 
20 

24 

27 
32 


Step  4  -  Calculation  of  One  v. a-        .  lation  Prora  VJholo     a.: 
to  a  '"as  Cube  at    (r  /'}• 

equation   (41 )  ]#    (Seo  Figure  X) 


rx/R 

r2A          f 

[» if 

*f^]^9A 

E.( 

^•(a 

Jl  rT.A     ->  / 

0.G0 

0 

3.5340 

0 

0.1 

3. 

o, 

677 

0.2 

.7798 

21, 

0.3 

4.1249 

» • 

^/  •  ■— 

1 . 7994 

-47 

0.3 

3.9762 

224 

O.G 

G.40 

458 

0.7 

,  D466 

k^j 

0.8 

3.3825 

61, 

0.9 

•  4312 

39, 

1.0 

1.904 

Then  by  numerical  integrations     j  fcoluian  4|    d(rg/fe)  -  51f188 


o^ 


-  '€■    - 


Stop  4  -  Continued 


TiA? 


o 

o.i 

0.2 
0.4 

o.r)r 


Tabulation  of   (q/A)    :as  to  cubos  at  all  radii. 


67,  . 

66,  sf?:: 

&  ,709 

62, 

59,085 

.55,691 

51,* 

45,210 

*  j.1»-> 

25,126 


71 

5890 
C704 
5473 

)9 
12 
35 

3470 

2743 

2425    (PARED) 

2215 


OJf- 


ep  5  -  Combine        'iation  Tarae 

(q/A)r  «    (q/A)sa+    (<JA)enlB8l<>a+  (qA)wall 


rx/h 

^3 

(q/A)  ,, 

uall 

(q/A)  ^  . 

v      amission 

(q/A)p 

0 

5971 

20 

-0794 

-2803 

0.1 

5090 

20 

-9557 

547 

0.2 

5704 

20 

-10107 

-  _ 

0.5 

5473 

21 

-10320 

-  335 

0.4 

5208 

22 

-10320 

-  )99 

0.5 

4909 

23 

-10127 

-5105 

0.6 

451S 

24 

-9751 

0.7 

5935 

25 

-9106 

-5095 

0.9 

3470 

27 

-3120 

-4631 

0.0 

2743 

32 

-6357 

-c 

0.05 

2425 

-4519 

-  074 

1.00 

2215 

-56 

+23 
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Fir 

us 
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XXT1 

r 

■  | 

.... 

im 

~ - 

it 

gdj! 

*j 

31 

'/ 

,* 

X 

% 

Vy 

* 

; 

n 

#• 

\ 

/ 

-— 

16 

— 

■e— 
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3£ 

/lr  ^ 

i* 

.... 

*»~ 

] 

' 

S* 

— 

— 



j 

f 

h 

~~ 

— ; — 

t 

- 

4fJ4 

—J 
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1 \ 

A»»o 

. 

-/-- 

•- 

— i 

; 

& 

j 

p 

--■ 

*••• 

•r   j 

| 

'   : 

- 

__ 

i 

\ 

! 

,  ..... 

1 
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— 

i 

- 

- 

.... 

'»»• 

/ 

•   • 

■   !■■■ 

...  .  . 

-„. 

-- 

--- 

.... 

— i— 

.. 

— 



; 
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— 

. 
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— fc 

-f— 

— 

* 



|i« 

X 

pi 

3 

*> 

m 

M 

lie 

<* 

:     . 

gj 

7 

** 
#. 

r 

£i 

£ 

* 

— 

M 

-v**o 

1 

f 

1 

% 

' 

*£ 

— 

— 

-- 

h 

- 

"ft 

! 
— 

— i — . 

•J«M 

' 

Mp 

— 

- 

_ 

i ' 

— i — 

1 — — 

— 

*4— i 

-I- 

-4— 

}— 

: 

— — 

i 



—- r 

' 

: 

^ 

j  ■■'■ 

i 

■ 
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i 

i 
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-  Plo\,  (q/A)  versuo  (y/r0) 
Draw  in  faro  eurv    e   '.  ~ 


snail  ccalc 
of      lal) 


(*V'  } 


(y/r0) 


c: 


(q/o. 


. 


• 


» 

0.    .'• 

0,70 

0.60 
0.50 

0.30 
0.20 
0.10 
0 


0 
C.3C22 

,  >527 

0.10 

o.i: 

v  .  20 

• 

OT  -\ 
»«J 

Q.C3 
0.." 

■0 
0.90 
1#00 


1#25 

2.50 

5 
10 
15 
20 
05 
2  . 
30 


5C9 


rl^ 
t1- 

t1'- 
t12 

*2 

-•  00 
-11 
-1'.- 

- 

_  ■ 

— Ol  . 

-  210 
- 

-4 
-609C 

-  •:. 
- 

-4035 

-  585 
-3C47 

-  BOS 


ps  7  and  8     -     obtain  the    (qA)     Integrals  in  t: 

tonparatur        Lstribu  aquations    (2 

(22)  an  ). 


^Tio  if    (0<yx.j5j_ 


\  (q/A)rd3T  ■  + 8,377     at  y+«  5  Bon»a 


Zono  II ,    (5<ys"<  30 )g   integration  by  planimotor 

+  •  y+ 

7  (q/A)r  r        (q/A)r  dy 


<*~7 


?r 


>  lG5.no" 

+  85.06 

-29.17 
-5^.00 
•64. 

-67.79 

-60. 7C 


Pr 


L_ 


10 

15 
20 

25 

20.45 
30 


0 


A  O^ 


^17G 

-5G2 
-466 


- 


laps  7  and         -     Continued 


Sono  III,      (50  <  y<    r   f,   Integration  by  planinoiier 


(y  /O 


*(   1  -  yVri   )/ 


o' 


<f^>* r  (o/a)„  dy* 


o  'J   W    yr(  1  -yVrp 


.  1527  -7  '.  0 

.10  -69. OS  -2 

.IS  -50.03  -8887 

-0.87  -63 

-4G.04  -6*; 

-42,69  -7 

0.40  -38.19  -10270 

0.50  -3l  .52  -12376 

0.G0  -37.34  -14437 

1.70  -40.61  -16G72 

0.30  -48.14  -101 

•0  -71.22  -224 

1.00  -CO  


Steps  9  anc1  10     -     Calculation  of  "Tov;  Temperature  Distribution 

equations    (20),    (22)  (24) 

idiation     Radiation  (q/A)0  T 

y       Integral     Coefficient  (3) &(£)  Terra        !    )-(&)    ( 
(Step  G)      (App,    IV) 


s  h 


® 

®      [ 
0 

(D 

© 

© 

0 

0 

0 

© 

0 

0 

0 

550 

3,  )( 

5 

♦8377 

0.  '0401 

-■ 

-34 

-"09 

959 

J. 0176 

10 

T 

>«  X3436 

♦14 

-se 

-r 

1217 

.  >352 

20 

+176 

+  6 

- 

-:.. 

150S 

5. 

.45 

-362 

- 

-7' 

-71 

1651 

27 

50 

- 

\ 

/ 

-16 

-749 

-733 

1672 

,10 

56.9 

•2006 

0.01718 

-34 

-160 

-120 

1700 

-  •  . 

-3827 

-66 

-260 

-194 

1366 

10 

-5377 

-92 

-• 

"Out/ 

-240 

1912 

- 

-6743 

-116 

-3. 

-271 

1945 

. 

-7980 

-137 

-4: 

— 1~, . 

1967 

.0 

-10270 

-176 

-505 

- 

2001 

0.5 

-12376 

-213 

-3  CO 

-347 

2019 

0.60 

-14467 

-249 

-606 

-557 

2029 

0.70 

•16672 

-236 

-044 

-; 

2030 

0.G0 

-19175 

-329 

-677 

- 

2020 

0.90 

-22447 

-386 

-706 

-  520 

1992 

1.00 

569 

V 

/ 



-733 



1950 
(FAI  ED) 

Ai'F 


CA 

cat  ion  of 

if 

r[       V//V 

r  ^ 

rx/R 

*1 

A. 

m-. 

.r. 
R 

a. — 

0 

1350 

24768 

0 

0.1 

1902 

20970 

30,52 

0.2 

2020 

28521 

G520 

. 

2030 

29090 

10104 

0.4 

2020 

29032 

14086 

0.5 

2010 

20405 

18118 

0.5 

2001 

27463 

21931 

0.7 

1967 

25645 

25507 

0.0 

1912 

22893 

27710 

0.9 

1790 

17903 

29003 

'  • 

55 

157 

03 

^->  ' 


The  Integral  of  column  4  with  respect  to  v^/u   from  0  to  1 

is  lGf153  by  Simpson's  Hulo 


(q/A) 


■■:: 


, 


(q/A) 


1 


(q/A)Wi     = 


1: 


4(x).(Ka!!)  s      (10,130)      3  1,024  Etu/nr.ft2 

?r 

-   E       »     -157     tu/hr.ft~ 


(q/A)w  ■+    (q/A)r 
S  € 


(q/A)        =  1GG7  Otu/ir.ft' 

Wp 


b  radiation  to  wall  from  the 
whole  "as 
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